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Abstract. In a previous paper (J. Phys. A 37 (2004) 9651-9668) we have given 
the Fuchsian linear differential equation satisfied by x^^\ the "three-particle" 
contribution to the susceptibility of the isotropic square lattice Ising model. This 
paper gives the details of the calculations (with some useful tricks and tools) 
allowing one to obtain long series in polynomial time. The method is based on 
series expansion in the variables that appear in the (n — l)-dimensional integrals 
representing the n-particle contribution to the isotropic square lattice Ising model 
susceptibility x- The integration rules are straightforward due to remarkable 
formulas we derived for these variables. We obtain without any numerical 
approximation x^^^ ^ ^ fuUy integrated series in the variable to = s/2/(l + s^), 
where s = sh{2K), with K = J/kT the conventional Ising model coupling 
constant. We also give some perspectives and comments on these results. 
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1. Introduction 



Since the formal expression for the magnetic susceptibility of square lattice Ising 
model, derived by T.T. Wu, B.M. McCoy C.A. Tracy and E. Barouch [T], and the 
exact expression of the Ising susceptibility written as an infinite sum: 

oo 

X{T) - (1) 

n=l 

of {n — l)-dimensional integrals [21 El El El El 0) the integration of the latter has 
become a long-standing problem in statistical physics. The above sum is restricted to 
odd (respectively even) n for the high (respectively low) temperature case. While x'^' 
is obtained directly without integration, and x^^^ is given in terms of elliptic integrals. 
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no closed forms for the higher terms are known. It is only recently that the differential 
equation for x^^^ has been found 

RecalUng the particle physics terminology of T. T. Wu et al., these (n — 1)- 
dimensional integrals are seen as successive "particle excitations", the first one, x^^\ 
appears to have been calculated by Syozi and Naya [5] for anisotropic Ising models (on 
the triangular lattice), and is actually a very good approximation to the susceptibility 
X, its anisotropic series expansion being also in remarkable agreement with the 
anisotropic series expansion for x- 

{l-ti){l-tl){l-tl) 

{1 — ti — t2 — t^ — tit2 — t2^3 ^ tits + ili2^3) 
16(1 + t lt2^3)(tl+^2t3)(t2 + ^lt3) (^3 + ^1*2) ^ 1/8 
{l-tmi-tl)Hl-tl) 

where ti, t2 and ^3 are the three standard high temperature variables, ti = th{Ki), 
for the triangular Ising model, and where M* is nothing but the magnetization of 
the triangular lattice Ising model for which the low temperature variables have been 
replaced by the high-temperature ones [TU| . 

For the square lattice Ising model, x^^'' is obtained as the ^3 = limit of JSJ for 
the anisotropic case, and as the ^3 = and t2 = ii limit for the isotropic case. The 
two-particle contribution x^"^^ is also known exactly and can be expressed in terms of 
the elliptic integrals E and K jj| (see also H82f) below) . These are the only two exact 
results on the n-particle contributions x^""*- 

Besides these two exact results on the n-particle contributions x^"\ the only two 
exact "global" results known for the whole susceptibility x^ correspond to the exact 
expression of the susceptibility satisfying the so-called disorder conditions and 
the exact functional relation corresponding to the so-called inversion relation |12| for 
the susceptibility HI)'. 

For one of the four disorder conditions of the checkerboard Ising model, ti + 
^2 ^3 ^4 — 0, the susceptibility of the checkerboard Ising model reduces to the simple 
rational expression 

X(il + t2 t, ^ 0) ^ {I^t2){l^t,)il+t,)il+t2t,t,) 

^ (X-t2tz){\-t^t^)[\-t2t^) ' 

or, more simply, on the anisotropic triangular lattice {t^ = 1 limitf): 

U , OA (l+^2)(l+t3)(l + ^2^3) ... 

Recall that Syozi and Naya algebraic expression Q for the anisotropic triangular Ising 
model reduces to the exact expression (0} when restricted to the disorder condition 
ti +^2^3 — 0. Therefore, one finds that all the x''"'''s {n > 1) of the anisotropic 
triangular! Ising model should vanish when restricted to the disorder condition. 

The so-called inversion relation 12^ for the susceptibility of the checkerboard Ising 
model and for the Ising model on the triangular lattice, read respectively |1L)I I13j : 

Xih, t2, ts, ti) + xi-ti, l/t2, -ts, l/ti) = (5) 
X(ii, ^2, t^) + xi-ti, 1/^2, -t3) = (6) 

f In this limit divide by a factor 2 to get the susceptibility per site. 

t This is also true for the four disorder conditions on the checkerboard Ising model, namely: 
x'"'(ti + t2tst4 = 0) = 0, n > 1. 
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to be combined with the obvious functional equations associated with the geometrical 
symmetries of the model (x(ii,i2,^3) — x{t2,ti,t3) — ■■ )■ It is straightforward to 
verify that the algebraic expression actually satisfies the previous inverse functional 
relation lO, and deduce, order by order in the n-particle expansion (Q, that all the 
n-particle contributions, x*-"-*, should also satisfy equations lO, together with the 
obvious geometrical symmetries. 

Of course all these ideas (x''^'' versus exact expressions of the susceptibility 
restricted to the disorder conditions, inversion relations together with geometrical 
functional equations, Kramers- Wannier duality, etc) can be worked out on the 
checkerboard Ising model (see for instance |10[ I14L This kind of analysis 

underlines the crucial role played by the modulus§ of the elliptic function 
parameterizing the checkerboard Ising model : 

p tUU+t,htl)il-tf) 2^/ 1 ^^^^ 

^.(i*+Wr)(l-i?) " \sh{2K,)-sh{2K,)J ^ > 

where the t*'s denote the dual variables, t* = {1 — tj)/{l + tj). For the anisotropic 
square Ising model, the inverse functional relation simply reads JOI (^3 = limit of 

my- 

xih, t2) + x(iAi, -^2) = (8) 

which can be combined with the obvious symmetry x(^ii ^2) = x(^2, ^i) • 

Using these last two exact identities, some re-summed high-temperature series 
expansions were performed for the susceptibility x of the anisotropic square 

lattice Ising model |1(JI I18L showing the occurrence of (2 fc -I- l)-th root of unity 
poles [ini, tl''^'^ = 1, for the {2k + 1) particle excitations x^^'' + ^^- This enabled 
A.J. Guttmann and I. G. Enting to conclude that x(iii ^2), as a function of t2, 
for ti fixed, has a natural boundary (the unit circle 1^2 1 = 1), and, thus, cannot 
be expressed as a solution of a linear differential equationjj of finite order N with 
polynomial coefficients of degree d, (even with N and d very large). Solutions of such 
linear differential equations with polynomial coefficients are called holonomic or D- 
finite |2Di 21 . Therefore, the susceptibility x is a "truly transcendental function" (in 
^2, ti fixed) if one takes for the definition of transcendental^ that it is not holonomic or 
D-finite. However, the x''"'''s, being multiple integrals of simple holonomic (algebraic!) 
expressions, they must be holonomic : they are actually solutions of linear differential 
equations. To sum up, x is not holonomic, but it is the sum of an infinite number of 
holonomic expressions. We have, here, a situation which is totally reminiscent of that 
encountered with Feynman graphs where all the individual Fcynman contributions are 
holonomic, but the sum of all these contributions is not holonomic (see for instance all 
the papers on multi-loop Fcynman integrals, nested sums and multiple poly logarithms, 
Euler-Zagicr sums, and other "motives" theory ,22, 23, 24, 25 ). In this respect, 
the analysis of the susceptibility x should be seen as the simplest example of such 
nested sum calculations, and, beyond the legitimate interest one has for the Ising 
model, should be seen as a "laboratory" for such calculations, hopefully displaying 

§ Recall that the spontaneous magnetization of the two-dimensional Ising model reads llbl 1171 : 
M = (1 - fc2)l/8. 

II For instance it cannot be expressed in terms of linear combination of special functions of 
mathematical physics like elliptic integrals or hypergeometric functions. 

^ For example the Gamma function is a transcendental function: it does not satisfy a linear 
differential equation, however it verifies finite-difference (linear and non-linear) functional equations. 
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the emergence of structures and symmetries. In this respect, x being considered as 
an infinite sum of holonomic n— particle excitations x*-"^ , a better knowledge of the 
analytical structure of the x^"^ 's is clearly a way to get a deeper understanding of the 
still mysterious analytic structure of x in the complex plane. 

Let us make here a few remarks. There are many "levels of complexity" among 
transcendental functions. Here we mean by "transcendental", non- holonomic, non-D- 
finite. As far as x is concerned, being transcendental, but expressible as the infinite 
sum of holonomic terms, one could seek solutions of simple, but non-linear, differential 
equations. For instance, it is worth recalling some particular§ non-linear differential 
equations, namely, the so-called Painleve trans cendentals |28[ I29| that are known to 
occur for some two-point correlation functions of the Ising models [SU] (Painleve 
transcendentals do not have natural boundaries). As a consequence of the previously 
mentioned natural boundary, the susceptibility x cannot be a solution of a linear 
differential equation, but it may well be a solution of a non-linear differential equation, 
even a very simple one! Actually, it is also worth recalling the Chazy equation |31[I32 |: 

—y- 2y—y-3(—yY (9) 
dw'^ dw-^ \dw J 

which arises in the study of third order ordinary differential equations having the 

"Painleve property" || PB^. This non- linear equation is probably the simplest example 

of an ordinary differential equation whose solutions have a natural boundary. For the 

Chazy equation (jSJ, the boundary is a movable circle% (its position depends on the 

initial conditions). The whole susceptibility x, being the sum of all the two-point 

correlation functions, one could, thus, perfectly imagine that Xj despite its unit circle 

natural boundary f , could well be a solution of a simple non-linear differential equation 

related to Painleve-like transcendants. 

While there are no algorithms to find such non-linear differential equations, 

implementing an algorithm to seek for a linear differential equation reduces to solving 

a set of linear equations. One knows that each of the x^"'''s is a solution of a 

linear differential equation, however, and unfortunately, the mathematical theorems 

that prove that the x*'"'''s are holonomic, because they are multiple integrals of 

holonomic expressions, do not give any upper bound on the order of the linear 

differential equation and no upper bound of the degree of the polynomials in front of the 

successive derivatives either [213 1^ • One can try to get enough coefficients of some 

high (respectively low) temperature series expansion to actually "guess" the linear 

differential equation satisfied by a given x''"'- Recalling the previous exact results 

that arise from a study of the anisotropic model, it might be tempting to generate 

anisotropic high temperature series expansions in order to use these exact results to 

get severe exact constraints on the coefficients of the series expansion. This can be 

done in principle, however, the number of coefficients in the two variables, and the 

combinatorial complexity associated with the finding of these coefficients, is too large 

§ They have fixed critical points I2til and are associated with some monodromy preserving property 
of some associated Schlesinger system 1271 . 

II The solutions have only poles for movable singularities S,§i. Movable means that the singularity's 
position varies as a function of the initial values. 

^ Heuristically, this can be simply understood introducing exponentially small terms in the expansion 
of the solutions (see (3.37) in 1341 '). or as Chazy did )31l 1321 . by relating its solutions to those of a 
linear hypergeometric equation, performing an infinite process of analytic continuations yielding a 
tesselation of the interior (respectively exterior) of a circle. 

f The susceptibility x of the isotropic Ising model has a unit circle natural boundary in the variable 
s = sh{2K), see JTOJ below and EIQ. 
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compared to the "advantage" one gets from the previous exact constraints. Therefore, 
we will ignore the previous exact anisotropic results, and will generate more standard 
isotropic series expansions for x^"^i ^-^d more specifically for x'^'- 

Recall that a huge amount of work had already been performed by B. Nickel [HUT] 
to generate isotropic series coeflticients for x''"' ■ More recently one should mention the 
generation by W. Orrick et al. |SS|, of coefficients of x using non-linear Painleve 
difference equations'^ for the correlation function EHl 133 EHl EH! ■ This second 
method, because it uses an exact non-linear difference equation was able to provide 
an algorithm for computing the successive coefficients in polynomial time |35) (namely 
0{N^)), instead of the exponential growth one expects at first sight. Roughly 
speaking, the first calculations by B. Nickel [HJ |7j were actually performed using 
the multiple integral form of the x*'"'''s and were thus able to provide series 
expansions for the x^^^'s separately. In contrast, the second method, which takes into 
account a fundamental non-linear symmetry, namely non-linear Painleve difference 
equations [371 03 EI j provides, as a consequence, a series expansion for the whole 
susceptibility X: where the various x^"'''s are difficult to disentangle. To sum up 
roughly, the first standard method is holonomic, or linear, oriented (decomposition in 
the holonomic x'^"'''s) while the second method is more non-holonomic, "non-linear 
oriented" (transcendental Painleve symmetries ...). 

As far as the singular behavior near the dominant ferromagnetic point T = 
Tc is concerned, it is numerically extremely hard to disentangle the |T — Td^/"*, 
\T-Tc\^/'^-\og\T-Tc\ and jT- T^p • log |T- contributions, and this was certainly 
a motivation to get longer series. This set of extensive results [HllZIESl (expansion up 
to order 116, then to order 257, and to order 323 for x) already enabled one to get 
a better understanding of the various singular behaviors. The short-distance terms 
were shown to have the form (T - TJp • (log \T - Td)'' with p > As far 

as the other singular points are concerned (non physical singularities in the complex 
plane), B. Nickel jH], in a remarkable paper showed that 

^(2n+l) jg g-j^ 

gular for the 

following values of s = sh(2J/kT) lying on the unit circle (k — m — excluded): 

2-{s + -) = + ^ + + 4r (10) 

y2n+i ^ -n < m, k < n 

the singularities being logarithmic branch points of order g2n(n+i)-i . ln(e) with 
e = \ — s/si where Si is one of the solutions of (|10|l . This confirms, this time 
for the isotropic model the existence of a natural boundary for x{s)^ namely \s\ = 1, 
and that x{^) is a transcendental function: it is not D-finite (holonomic) as a function 
of s. 

Our approach here is clearly of holonomic type, underlining the crucial role played 
by some well-suited hypergeometric functions. The emergence of these structures 
actually allows us to get extremely efficient polynomial time calculations (namely N"^^ 
for x^'^^)j which enable us, in particular, to calculate the series coefficients separately ^ 
and not from a recursion that requires the storage of all the previous data. In 
the sequel, we give a fully integrated expansion of the three-particle susceptibility 
contribution x^'^^ as multisums of hypergeometric functions. This expression is used 
to generate a long series in the variable w = s/2/(l + s^). With this long series, 
we succeeded in obtaining the homogeneous seventh order Fuchsian differential 
equation satisfied by x^^^- 

tl More generically this corresponds to the notion of the so-called Hirota bilinear equations. 
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This paper is organized as follows: in Section 2, we present the basic features of 
our expansion method that allow us to obtain the fully integrated x'"^^ ^ multisum 
of products of three hypergeometric functions, without any numerical approximation. 
The three-particle susceptibility x^^^ is written, for series generation purposes, as 
the sum of a closed, but involved, expression (sum of linear, quadratic and cubic 
products of elliptic integrals), of a second simple (non closed, for the moment, but easy 
enough to compute) sum, and a last sum which is much harder to calculate, and which 
requires most of the computing time. The variable w being well-suited to both high 
and low temperature regimes, the expansion method we use is generalizable to other 
X^"-''s {n odd, or even). In Section 3, we recall the homogeneous linear differential 
equation and discuss the new singularities discovered in [S]. In Section 4, using a 
differential approximant method, we show how the new non apparent singularities 
can be discovered even when the series are not long enough to obtain the exact linear 
differential equation. Finally, Section 5 contains our conclusions. 



2. The fully integrated expansion of x^^^ 
2.1. The expansion method 

Let us focus on the third contribution to the susceptibility x defined by the double 
integral as given in 6 : 

X^^\^) = {l-^'Y'Vs-x'^''\s) 

x^'Hs) - ^ r dc^, rdc^, • yiy.y3 • (^^li^i) . ^(3) (u) 

47r^ Jo Jo \l-x1X2X3J 

with 

= \ — ; : n. 9 ' (12) 

1 + — s COS (pj + ■^/(l + s"^ — s cos (f>jY — 
Vi, = r— „ , ,„ „ , i = 1,2,3, 01 02 + (/'a = 0. 

-I- — SCOS(f)jY — 

Many forms for H'^'^^ may be taken [HI [7] and are equivalent for integration purposes, 
e-g-, 

i?^'^ = /23-(/3i + 4^), = (sin0. - sin0,) - (13) 

Instead of the variable s, we found it more suitable to use w = ^s/{l + s^) where 
the unit circle \s\ = 1 becomes ] — 00, —1/4] U [1/4, oo[ in w. This variable allows to 
deal with both limits (high and low s) on an equal footing. 

From now on, we will use the scaled variables: 

a; . :^ 51 :^ ^ (14) 

^ w 1 - 2w cos (j>j + ~ 2wcos(j)jY - 4w^' 
1 



2w ^(1 - 2wcos(/)j)2 - 4^2 
which behave like: 1 -|- 0{w) at small w. These variables are related by: 



Vj • (1-w x^) = X,. (15) 

Performing the integrals in Hll|) is highly non trivial. For this, one may expand 
symbolically the integrand in ()11|) in the variable w, and integrate the angular part. 
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This way one faces some 18 sums to carry out. These sums come from the expansion 
of the quantities Xj and yj, and also of the denominators. The convoluted character 
of the integrand, and especially 1/(1 — w^xiX2X3), makes the symbolic expansion in 
w almost intractable. 

To overcome this difficulty we expand the integrand of x^^"^ i^i the various variables 
Xj . The expansion depends now only on combinations of the form : 

yiXi' ■ V2X2'' ■ VS^s^ ■ -P«i,«2,n3(<^l>2, 03) (16) 

where Pni,n2,n3{4'ii 4>2i4>3) is a polynomial in circular functions of the angles <j)j. The 
structure (|15|l appearing in the integrand, after the expansion in the variables Xj, is 
independent of the particular form taken for H^^^ . 

We succeeded in deriving a remarkable formula for UjX^ that carries only one 
summation index and reads (we drop the indices) : 

00 

a{k, n) cos kip (17) 

fe=i 

where a{k,n) is a non terminating hypergeometric series that reads: 

aik,n)=('^y (18) 



4^3 



f {l+m) (1 + m) (2 + m) (2 + m) ^ , ,^ ^ , 2\ 
[ ^ , ^ , ^ , ^ ;l + fc,l + n, l + m;16w j 

where m ^ k + n. Note that a(fc, n) ~ a{n, k). 

With this Fourier series, the angular integration of the form (|16f) becomes 
straightforwardi . 

The method described above being independent of the particular choice of H'^^\ 
let us consider the form Noting that : 

XiXj _ ^ Xi — Xj 



h = -77{x^- X,) COt(0,/2 + (f>j/2). (20) 



1 — XiXj 2 COS 0i — COS 

the fij 's can be rewritten as : 

w 
2 

H^^^ becomes simply a polynomial in the x^'s : 

= — . (21) 

-{X2 - XsY ■ COt^(y ) - {Xl - X3){X2 - X3) ■ COt(y) • C0t(y ) J . 

Expanding the integrand of x*-"^-* in the variables Xj : 

^(3)(^)^£^/ dc^^ I dct,2- (22) 

00 

E 

n=l 

X Note that we can deduce from the Fourier series 1171 that yx" also has an expansion in the variable 
w (respectively siT?((f>/2)) with coefficients as hypergeometric functions of the variable svc?(4>/2) 
(respectively w). Actually, the expansion I17i is more suitable for angular integration purposes. 



47r2 

2/12/22/3 • (1 + 2 ^ w^"(a;iX2X3)" 
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together with the expression (EU for H'-^\ one notes that the integrand depends only 
on combinations of the form : 

-yax^^ -yaa^a' •cot2('^i/2) (23) 
yix^ • y2xT ■ ysa;^ • cot((/)i/2) cot(02/2) (24) 

which have straightforward integration rules. 

The problem of integration is thus settled with a limited number of sums. In 
Appendix A we give the integration rules used and show how the cancellation of the 
"artificial" singularity introduced by taking the particular form ^21^ for i?^^^ occurs. 

The expansion of the integrand in the variables Xj , together with the remarkable 
formula Ijl?!) . allow us to obtain x*^'^' (w) as a fully integrated expansion (see Appendices 
B and C), namely: 

X^^^w) = 8 w;^ • (xif (u-) + 2w^- SH) (25) 
where are sums given in Appendix B, and Xd i'^) ^ closed expression: 

*^'^(^)= 256(1 - 4^/^3- (26) 

3 i 2 i 

1=0 j=Q i=0 j=0 

K and E being functions of w related to elliptic integrals K and E: 
k = 2^1(1/2, 1/2; 1; 16^2^ 2K{Aw)/t:, 

E = 2Fif-l/2,l/2;l;16u;2) = 2E{Aw)/t:, (27) 



and where the Qij^s and Pij 's are polynomials in w given in Appendix D. One should 
note that this closed expression H26I) is not " natural" : the separation in (|25|1 is not 
unique. It comes from our ability to evaluate some sums, and is only made for series 
generation purposes (see Appendices B and C for details). In fact, x^^) ||25|) can be 
written as a multisum of product of three hypergeometric functions a(fc, n), as can be 
clearly seen§ in p6|) . 

One recalls that the key ingredient in obtaining H25|l is the Fourier expansion 
of ya;", a quantity appearing in any x*-"''. The method is thus generalizable, 
independent of whether n is even or odd, since the variable w has, by construction, 
Kramers- Wannier duality invariance (s <~-> 1/s). Appendices E and F show how the 
known result ITI for x^^^ appears, in this framework, as a sum of products of two 
hypergeometric functions merging into a single hypergeometric function. 

For the x^'^^ case, we see that the mechanism oi fusion relations on hypergeometric 
functions that we saw previously for x^'^\ may also hold similarly. However, we have 
not discovered a similar fusion mechanism for x'-^'' • 



2.2. Series generation 

From our integrated forms of x^^"*: the generation of series coefficients becomes 
straightforward. Recall that (|25ll is already integrated, and, thus, the computing 

§ For this multisum of product of three hypergeometric functions can be seen in the definition 

of the auxiliary functions (see Appendices A, B and C) in terms of a(fc,n). 
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time to obtain the series coefficients comes from the evaluation of the sums. For the 
forms used, this time§ is of order N'^. 

We have been able to generate a long series of coefficients from the expression 
H25|l up to order 490. More precisely, the series expansion for S in (|25|l . was, in fact, 
obtained as the sum of two series, a "hard to compute" series up to order 490, S/j, 
requiring most of the computing time||, and a simpler one, S^, requiring much less 

(3) 

computing time, but of course, more time than the closed expression Xcj • They read 
respectively : 

2w* ■ Ehiw) = 16w'^ + 88w^ + lOOSw^ + 5144w^ + • • • (28) 
2w^ ■ Esiw) =^ iw^ ~ 4w^ + 6w^ ~ 1628u;^ - 11738u;^ + • • • (29) 

and ^ : 

x''^\w) = 1 + 36w2 + Aw^ + 864^4 + n2w^ + 17518u;'5 + • • • (30) 
The expansion of x'-^'^ thus reads : 



^■^^^ = 1 + 36^2 + + SSAw"^ + 196w^ + 18532u;^ + • • • 
4309904156072449034005168913683182133673083305758179058 



(31) 



1338044228506645721942836227131943916427611183965848882 

0854328108563517776050668285313206447166513353732686903 

6438360433784746161733072061779435657865203525518535003 

4549334011001139634569887672761933044541359314192379119 

5381220947969620- w^^^ + ••• 

in agreement with the coefficients obtained by B. Nickel and those obtained by 
A. J. Guttmann et aLf , these last expansions being, in fact, written in terms of the 
u = s/2 variable : 

^ — )—L = 1 + 4^3 + 16u^ + 4u^ + 20u^ + 84u^ + 247u** + • • • (32) 
4^8 

For large values of TV, the coefficients Cat, on the right hand side of H31() and H32() . 
grow, in the variable w as: 

Cjv~ 13.5x4^, iVeven, ~ 11 x 4^, iVodd (33) 

and in the variable u = s/2, as C^-* ~ 15 x 2^. 

Our long series is obtained as the sum of three contributions of different 
algorithmic complexity. The contribution H3UI) being given by an exact closed form, 
we were tempted to seek a linear differential equation satisfied by x^^ . Appendix 
H contains a report on this matter, which suggests that seeking for the differential 
equation for x'"^\ seen as the sum (|25|) . is in fact much simpler than seeking the 
differential equation for each constituent. The whole sum (|25|l is, thus, a better 
candidate for seeking for a linear differential equation. 

§ The overall computing time can be reduced by using the inhomogeneous recurrences on, and 
relations between the coefficients of the functions introduced (see Appendix G). 

II Expression Sj, actually corresponds to the true limitation in our series expansion. It corresponds 
to the third term at the right hand side of 1631 . see Appendix B below. 

1 This sum is actually analytic near to = 0, while the expansions of the various cubic sums, quadratic 
sums and linear sums in K and E in the closed form 1261 exhibit l/w^^, ■ ■ ■ poles, 

t Private communication : W. Orrick and A. J. Guttmann obtained the series up to 250 coefficients, 
in agreement with our 490 coefficients. 
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3. The Fuchsian difFerential equation satisfied by x'^^ 

With our long series, and with a dedicated program, we have succeeded in obtaining 
the differential equation for x^^^ that is given in |S] and recalled here: 

E «»(^) • ^^(«^) = (34) 

n=0 

with : 

aj = (1 - w) (1 + 2w) (1 - 4w)^(l + iwf {l + 3w + W) P7{w), 

ae ^w^{l-4w)^{l + AwfPe{w), ag ^ w^l - 4wf {I + Aw) P^iw), 

= w4 (1 _ 4u,)2 p^{w), 03 = w^{\- Aw) P^iw), 

a2^w'^P2{w), ai=wPi{w), gq ^ Pq{w). (35) 

where Pji^w), Pq{w) ■ ■ •, Poiw) are polynomials of degree respectively 28, 34, 36, 38, 
39, 40, 40 and 36 given in pi. 

The differential equation H34|) is an equation of the Fuchsian type since there are no 
singular points, finite or infinite, other than regular singular points. The singularities 
correspond to the roots of the polynomial in front of the highest derivative in 134|l . 
The ferromagnetic (w = 1/4 or s = 1), antiferromagnetic {w = —1/4 or s = — 1), 
zero or infinite temperature (w — or s = 0, oo), non-physical (w = ooors = ±i) 
and Nickel singularities [w — 1, —1/2 or 2 — s + 2s^ = 0, 1 + s + = 0) are indeed 
regular singular points of the differential equation. The last two are given by Nickel's 
relation pO|l for n = 1 . Besides the known singularities mentioned above, we note 
the occurrence of the roots of the polynomial Pi(w) of degree 28 in w, and of the two 
quadratic numbers l-|-3w + 4u;^ = which are not of Nickel's form H10|l . 

Applying Frobenius's method [211, we give in Table 1, for each regular singular 
point, the critical exponents, the number of solutions in Frobenius form, and the 
maximum power of the logarithmic terms. One sees that only log and log^ terms 
appear in the local solutions of 134(1 . 

Table 1: Critical exponents for each regular singular point. Np and P 
are, respectively, the number of solutions in Frobenius form, and the maximum 
power of the logarithmic terms for each singularity of H34|) . is any of the 28 
roots of Pj{w). We have also shown the corresponding roots in the s variable. 
In the variable s the local exponents for w = ±1/4 are twice those given. 
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One can see that the roots of Py are apparent singularities, the local solutions carry 
no logarithmic terms and are analytic since the exponents are all positive integers. 
Consider the apparent singularity effect from the series expansion x^^"^ viewpoint 
of this paper. If the roots of Py were not apparent and were singularities of x^^\ 
the coefhcients of (|31ll . instead of being dominated by the w — ±1/4 singularities, 
yielding a 4^ growth of the coefficients, would be dominated by the smallest (in 
modulus) root of Pj, namely w ~ —0.0424, yielding a growth ~ (23.553)^ of the 
coefficients. We have checked this on many (slight) deformations of (|34|l . keeping 
all but one of the P,i's unchanged (in particular Py), and slightly modifying one P„ 
(n ^ 7). One immediately sees a change in the growth of the coefficients, with the 4^ 
growth changing drastically into a ~ (23.553)^ growth. It is thus crucial that the roots 
of Py are apparent singularities, which require some very specific "alchemical tuning" 
between the various P„'s. In other words, the Fuchsian equation H34|l corresponds to 
a rather rigid structure, since the apparent singularity character of Py is an unstable 
property. 

The two quadratic numbers, w = (—3 ± i\/7)/8 (or 1 + 3w + 4w^ = 0), which 
read in the s variable 

(2s2 + s + l)(s2 + s + 2) = (36) 

have their roots lying respectively on the circles |s| = l/\/2 and \s\ — \/2. Near these 
points, one of the solutions of (|34|) has weak singular behavior (a log singularity). 
Note however, that some of these logarithmic terms, occurring in the general solution 
of H34|l and shown in Table 1 may not exist^f^ in the particular solution x^^^ . Actually, 
the two unexpected quadratic numbers correspond to singularities of solutions of (|34|l 
which behave locally like (1 + 3w + 4 w^) • ln(l + 310 + Aw'^) ■ qi + q2, where qi and 

^ See, for instance, the two simple rational and algebraic solutions Si and S2 below (see 1371 ) which 
are free of the 1 — w = 0, l+2io = 0, w = and 1 + 3w +4 w'^ = singularities. 
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q2 are analytic functions near the two quadratic roots 1 + 3w + Aw"^ — 0. This weak 
singularityf gives a 2^ contribution to the growth of the w-coefRcients. However this 
does not mean that each solution| of H34|l has such {l + 3w +4 w'^) • ln(l + 3w + 4 w^) 
singular behavior near the two roots 1 + + 4 = 0. Actually, as far as the 
"physical" solution x^^^ is concerned, and as far as the series analysis view-point 
developed in this paper is concerned, this (1 + 3w + 4 w^) • ln(l + 3w +4 nP') singular 
behavior should be excluded : in the s variable, instead of the w variable, it would give 
(beyond other terms) a (1 + s + 2 s^) • ln(l + s + 2 s^) term, yielding a a/2 growth, 
which would be in contradiction with the growth of the coefficients of x!"^^ ™ the s 
variable. This requires further investigations that we will address in a forthcoming 
publication. 

To end this section and to be complete, we will recall the pertinent algebraic 
properties of the Fuchsian differential equation whose analysis was sketched in jSj. 
The Fuchsian differential equation (|34() has two remarkable rational and algebraic 
solutions, namely: 

Si{w) = S2{w) = (37) 

^ ' l-4w' ' (1 -4u;) VI- 16U.2 ^ ' 

This results in very important factorization properties for L7, the seventh order 
linear differential operator, corresponding to the Fuchsian differential equation H34|l : 

L - — + — Ya.— (38) 
dw'^ aj ^ dw^ 

' k=0 

The adjoint differential operator L7 has the following rational solution§: 

o*f ^ /H ■ Qejw) , , 

" — 3 p / N (39) 

where /(uj) ~ (1 — u))(l + 2u')(l — 4?ij)^(l + 4u))^(l + 3ti; + 4?«2) and Qg is a polynomial 
of degree 28 given in 

All these findings imply the following factorization of L7 : 

L7 = Me -Li, L7 ^ Ne-Ni (40) 

L7 - Ml • Le, Le = L^- iVi (41) 

where Mg, Nq and L5 are operators of order respectively six, six and five, and where 
the order-one operators Li, A^i and Mi are such Li(S'i) = 0, A^i(5'2) = and 
Ml (5*1) = 0. The operators Lq, Mq, Nq and L5 can easily be calculated from the 
previous factorization relations, by right or left division by the first order operators 
Li, Ni and Mi defined above. The decompositions (|40l 14 111 imply that the general 
solution of (|34() is a linear combination of and a solution, S{Lq), of the sixth 
order linear homogeneous differential equation associated with Lq. The particular 
"physical" solution of L7, x^^\ can be written as x^^^ = a ■ Si + S{Lq), where a is 
deduced from Lq(x'^^"' — a- Si) = and has the remarkably simple value a = 1/3. The 
homogeneous operator Lq has the polynomial f{w) Qaiw) in front of the highest 
derivative. The roots of the polynomial Qq{w) are also apparent singularities of the 
linear homogeneous differential equation associated with Lg. 

t This "weakly" singular behavior can also be seen in the monodromy matrix (see 8 ) associated 

with these two roots, where one finds a nilpotent matrix of order two (no log^ term). 

X A solution of 1341 can exhibit only a restricted set of the whole set of singularities of the Fuchsian 

equation: this is very clear with Si, or S2 (see I37i below). 

§ We thank Jacques- Arthur Weil for the remarkable results 1391 and 1411 . 
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4. Towards the new non apparent singularities: a Diff-Pade method 

From the ideas developped in Section 2, we were able to obtain a very long series 
expansion for x*-'^-', sufRciently long to actually identify the Fuchsian equation for x^'^''- 
However, when considering further similar calculations on the next x''"^'s (x^*'') X^^\ 
...), we may not be able to obtain large enough series to identify the corresponding 
Fuchsian equations. One may ask the following question: what kind of results can one 
get from the analysis of long series, but not long enough to find the exact Fuchsian 
equation. This obviously supposes taking into account some of the ideas previously 
encountered, namely the homogeneous character of the linear differential equation, the 
Fuchsian character of some known singularities (like s = ±i and Nickel's singularities, 
see ()10|l ). the possible occurrence of unexpected, but simple, regular singularities (like 
the quadratic numbers l+3i(j+4w^ = 0), and above all, the possible occurrence of 
quite a large set of apparent singularities\\ (like the 28 roots of polynomial Py). 

In fact the preliminary studies we performed before finding the Fuchsian equation 
give some hint as to the kind of results one might expect ^. 

One knows that x^'^^ satisfies a linear differential equation. For a given number 
of terms N in the x^^^ /w^ series, there is a linear differential equation of order q that 
reproduces the first N — q terms but may fail for subsequent coefficients. Let us show 
what kind of information arises from the differential equations obtained for increasing 
order g, and increasing degrees of the polynomials in front of the derivatives. We will 
consider the following form : 

n=Q 

There are various ways to choose the degrees of the polynomials p„(w)'s. One may 
take: 

degipn) = deg{pn-i) + 1, deg{pn) = n + n. (43) 

The differential equation (|42|l then has l/2(g+ l)(g + 2/i + 2) unknowns. The number 
of terms in the x*-'^-' /w^ series used is: 

N = i(q + l)(q + 2/i + 2)-l-g. (44) 

One should note that the results below are not dependent on the form H43|) we took, 
and are also obtained when the degrees of the polynomials Pn{w) are taken to be 
equalf . 

We systematically considered linear homogeneous differential equations with 
polynomial coefficients (in w) of increasing orders. The degrees of the various 
polynomial coefficients are taken as large as possible. With such a strategy the order 
of the homogeneous linear differential equation, and the degrees of the polynomials, 
are not on the same footing (the order is much smaller than the degrees and the order 
of a differential equation is a much more fundamental character of the differential 
equation). For N terms in the series, one then explores all the linear differential 
equations reproducing these terms with {q,fJ-) staying on the hyperbola H44|) . 

II Which suggests considering Fuchsian equations with quite large degrees for the polynomial 
coefficients, compared to the order of the homogeneous linear differential equation. 
1 At this step one can certainely consider using the differential-Pade approximation calculations 
developped by A.J. Guttmann et al. 1401 . or similar differential-Pade approximations considered by 
other authors 1411 . 

t The form 1431 assures -u; = oo as a regular singularity of the differential equation. 



Series expansion method and ODE of x^^^ 



14 



Let us fix g = 4 and compute the linear differential equation for increasing values 
of /i. The polynomial in front of the highest derivative is solved for the variable w. 
For /i = 6, one gets w = —1/2 with three correct digits, w = — 1/4 (double) with five 
digits, w — +1/4 (double) with four digits, and w = 1 with two digits. With these 
values of w, we get four other solutions. As fi increases, the number of correct digits of 
the roots w = —1/2, w = —1/4, w — 1/4 and w = 1 increases. For /i = 15, the values 
become correct to, respectively, 8, 13, 14 and 6 digits. For ^ = 28, they are correct to 
14, 20, 26 and 13 digits. The roots w = ±1/4 still have multiplicity 2. One observes 
that for q — 4: and n = 28, the number of series coefficients used is just 150 terms. 
Further, as fj, increases, the other solutions increase in number {pq(w) has q + fJ- roots) 
and do not converge to stable values. At around fi ^ 19 a, root oscillating around 
— 1/4 appears. One can imagine that the root w — —1/4 is going to be of multiplicity 
3. For the root w — +1/4, one has to wait until ^ = 28 to see this happen. 

The differential equations we obtain are quite cumbersome and carry as much 
information (in the number of unknowns) as the TV terms used. This is, then, just 
another way to encode the information contained in the N terms. What is remarkable 
in this procedure is the appearance of the singularities with increasingly correct digits 
and the fact that they converge to the correct multiplicity as ^ increases. What is more 
remarkable is the appearance of the other unexpected singularities 1 + 3w + iw^ = 
which were obtained with three digits at {q = 4, /i = 2), and with nine correct digits 
at (g = 6, /i = 35). 

Within this systematic approximation scheme, we actually discovered, at smaller 
orders than the actual one (namely order seven), the existence of the unexpected 
quadratic numbers singularities 1 + 3w + Aw^ = 0, with an extremely good level of 
confidence. This indicates that such a systematic procedure was probably converging 
towards some exact result, which we actually found with the next orders. This method 
is obviously applicable to any series that satisfies a linear homogeneous differential 
equation. Actually we will show, in a forthcoming publication, that we can actually 
find the singularities of the successive even when we are not able to find the 

ODE'S satisfied by the x^"^'s. 

5. Conclusion 

Considering the isotropic Ising square lattice model susceptibility, we focused on the 
third order component x^^^ in the form given by Nickel [S]. We used an expansion in 
the variables xj that appear in the integrand of the double integral defining x*-^^ instead 
of the variable w (or s). The angular integration becomes straightforward thanks to 
the remarkable formula we derived for the quantity yj . This formula is a Fourier 
series carrying one summation index and where the coefficients are non terminating 
4F3 hypergeometric functions. We succeeded, in this way, to write x^^^ as a fully 
integrated expansion, containing few sums of products of three such hypergeometric 
functions. For x^^\ the sums bear on products of two hypergeometric functions but 
due to remarkable identities, they "fuse" into a single hypergeometric function giving 
the known result (see Appendix E and F). We note that these fusion identities on 
hypergeometric functions occurring for x*^^-* , may be the sign of deep structures and 
symmetries. We have not yet recognized a similar mechanism for x^^\ but it seems 
plausible. 

We used the fully integrated multisum giving x''^'' to generate a long series 
of coefficients with a polynomial time {N'^) algorithm. We recall that the series 
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coefficients are obtained without any numerical approximation. We were able to obtain 
the series up to order 490 in the variable w, and used it to find the seventh order 
linear differential equation (|34|l satisfied by x'"^^ ■ This limitation to 490 terms is a 
reflection of the minimal computer resources used. Much longer series can readily be 
obtained. This differential equation, which is of Fuchsian type, is highly non trivial 
and structured. This result shows that quite involved sums of products of three 
hypergeometric functions can be solutions of quite simple linear differential equations. 

The ferromagnetic, antiferromagnetic, zero or infinite temperature points, the 
non-physical and Nickel singularities are all regular singular points of the differential 
equation. The Fuchsian differential equation shows other simple singularities, which 
are not of Nickel's form, and are lying on the circles \s\ = \/2 and \s\ ~ 
We also note the occurrence of a polynomial of degree 28 in w whose roots are 
apparent singularities of the Fuchsian differential equation. We also saw that any slight 
deformation of the Fuchsian differential equation drastically modifies the growth of 
the coefficients in the series expansions of the solutions, confirming the highly selective 
character of this ordinary differential equation. 

The behavior of the local solutions of the Fuchsian differential equation at each 
regular singular point is sketched and shows the occurrence of logarithmic terms up 
to log^. The question of whether these terms occur in the physical solution x^'^^ is of 
the utmost importance and, particularly, for the new singularities 1 + 3w + Aw^ = 
not previously known|. We will report on the integration constants corresponding to 
the physical solution x^^^ in a forthcoming publication. 

The Fuchsian differential equation has two simple rational, and algebraic, 
solutions enabling one to give the factorization of the differential operator 
corresponding to the Fuchsian differential equation. The solution is a linear 
combination of a rational solution, namely w/{\ — 4^), and of a solution of a sixth 
order homogeneous linear differential equation. 

Recall that the variable w we used deals with both high and low temperature 
cases on an equal footing. Furthermore, the key ingredient in obtaining (|25|) . and 
finally H34() . is the Fourier expansion oi yx", a quantity appearing in any x^"^- So 
the expansion method we used is not specific to the third contribution x^'^^ can 
be generalized, mutatis mutandis, to the other x''"\ and in particular the ones with 
n even, associated with the low-temperature susceptibility (in Appendix E we discuss 
the form of x*-^-* in this variable). The next step is clearly the evaluation of x'-'''' which 
can be seen [HIE] to be a function of s'^, and, in fact, a function of w'^ . This gives 
some hope of finding the corresponding Fuchsian differential equation. 

The Fuchsian differential equation, satisfied by x^'^\ is certainly an important 
step towards the understanding of the analytical structure of isotropic Ising square 
lattice model susceptibility. Along these lines, one can easily imagine finding the 
Fuchsian differential equation satisfied by x'^^'' for the isotropic triangular Ising model. 
However, recalling the nice and simple exact results we displayed in the introduction 
for the anisotropic square Ising model, on the anisotropic triangular Ising model, and 
on the checkerboard Ising model, it is tempting to seek for the homogeneous "Fuchsian 
partial differential equation" (generalizing our Fuchsian differential equation) in the 
corresponding two, three or four, high temperature variables ti, t2, and ^4. No 
doubt the modulus k of the elliptic functions parameterizing the Ising model (see 
((TJ) should play a special role in writing down a simple enough expression for this 

t For 1 + 3w) + 4^2 = singularities, in Section 3, we gave arguments in favor of its absence in x'^'- 
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Fuchsian PDE. Having such a Fuchsian PDF, it will be extremely interesting to 
see how this homogeneous Fuchsian PDF can be invariant§ (or covariant since it 
is homogeneous) under the inversion relation symmetries ® , ||SJ| and the whole set of 
birational transformations j43ll44lE5j generated by these inversion relation symmetries 
combined with the geometrical symmetries. 

Acknowledgments We thank Jacques- Arthur Weil for valuable comments and 
calculations on our Fuchsian equation. J-M.M would like to thank A. J. Guttmann 
and W. Orrick for extensive discussions and for a large number of e-mail exchanges 
and calculations more recently. We would like to particularly thank A. J. Guttmann 
for careful checks on our scries. (S. B) and (S. H) acknowledge partial support from 
PNR3- Algeria. 



6. Appendix A 

In this Appendix, we give the various formulas used to obtain H25() . We use < • • • > to 
mean the normalized angular integrations Yii "^^^ variables x and y are defined 
as in (|14|l . From (|17|l . one gets: 

< yx" coa{p(t)) > ~ a{n,p) — w'' a{p,n). (45) 

Recall that a{n,p) is defined in (|18|) in terms of 4_F'3 hypergeometric functions. 
Define the following functions: 

gin^p) = V(i + 1) • w'a{n + i + l,p) = w^p < — cos(p0) >, 

^ {l-wxy 

v(n,p) — w^a(n + i + l,p) = w^^ < cos(p0) > . (46) 

^-^ 1 — wx 

1=0 

From the integral representations, one has: 

a{n + l,p) ^ v{n,p) ~ w v{n + l,p), 
v{n,p) = g{n,p) -wg{n+l,p). 
The following integrals can, then, be computed 

< yx""- cot(0/2) sin(fc(/)) > yQX^ - w'" ■ (^a{k, n) + 2wv{k, nfj , 

< ya:" cot((/)/2)^ cos{k(f)) > = - 2fc yQX^ - w'' (^a{k, n) + Awg{k, n) 

1 

sin(0'/2) 

where xq and j/o a-re the variables defined in H14|l and taken at = 0. 
Some integrals appearing in the intermediate steps are {k > 1): 

l~cos{k(j)) ( 2(k-p), k>p 

< sin(0/2)^ ^ =lo' k< p (48) 

1, k > p 

< cot((/)/2) sm{k(j)) cos(p0) > ^ { 1/2, k = p (49) 



k <p 



§ See for instance 1421 . 
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We now turn to the "artificial" singularity appearing on the right hand side of 
H47|l . This term contributes to the total integrand Hll|l as: 

< ^ ("^^^^ + ^f:^'' dip)) with (50) 



p=l 



dip) =yoxP(^Sip + 2,p)-Sip+l,p+l)), 

OO 

S'(p, g) = a(p, 0)a(g, 0) + 2 ^ ly^^ a(p, A:)a(q, fc). (51) 



fe=i 



Noting that 

Sip, q) = < y^x^+'i > = <yxP yx" > (52) 

it is straightforward to see the cancellation of dip) and, thus, of (|5U|) as it should. 

Note that formula (|17|l has been used extensively in the following appendices, and 
each time we have an expression similar to (I52II , we manage to write it in terms of H17(l . 
Another trick used throughout the appendices is the technique of switching to-and-fro 
between the sums and the integral representation, and the judicious selection of the 
correct place to use ((T^ . 

7. Appendix B 

Expansion of Hll() in terms of yja;", using H17|l and the integration rules given in 
Appendix A, gives for x^^^^: 

oo 

X^^'> = 8 (/(0) + 2^u;3f/(p)), (53) 



with 



Hp) = a < y^y^y^ ixiX2X3)p h^^^ > 

oo 

= Cip) + w^FiO,p) + 2w^ J2 w^^Fik,p), (54) 

k=l 

where H'^^^ is taken in the form 121|) . Fik,p) consists of sums of cubic terms of the 
hypergeometric function (|18|l . and can be written in terms of this hypergeometric 
function and the auxiliary functions gik,p) and w(fc,p) defined in Appendix A, as: 

Fik,p) ^gik+l,p) (^aik,p + 1)^ - aik,p + 2)aik,p)) 

+ a(fc,p + 2) vik,p)'^ + aik,p) vik,p+ 1)^ 
-2aik,p+l)vik,p)vik,p+l). (55) 
The function C(p), written as: 

C(P)^ -I yl4^ ip) + y^xl C2 ip) , (56) 
consists of a linear term in the a's: 

Clip) = a(0,p + 2) + a;2a(0,p)-2xoa(0,p+l). 
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and a quadratic term in the a's: 

C2(p) = (To(p) + 2u;cri(p) + 2wa;o (T2(p), (57) 
where the sums (Ji(p) (i = 0, 1, 2) are given by: 

oo 

'^ob)- ^(fc + l)u;'^-(a(fc,p+l)' -a{k,p + 2)a{k,p)^, 

k=l 

oo 

CTi(p) = ^w^'' {a{k,p+l)v{k,p+l) - a[k,p + 2)v{k,p) 
fc=i 

oo 



fc=i 



To get the final form of x*-^-* , further manipulations are made on that part of 7(0) 
and I{p) in (|53|l containing the C"s. The aim is to get only one summation instead of 
two. We show in Appendix C how to obtain the part of ^b'6\ containing the C's. It 
reads: 

oo 

C(0) +2^w3PC(p) = Ci+c^c where (58) 

3 

6 = yo C(0) + ^ • (jo(0) - a(l, 0)) (59) 

oo 

ac = 2 2/3 ^ w;3«-2^„-i - u; (1 - u;a;o) • Ji(n)) 



n=l 

where >/o(?^) = a(0, ti)^ — w^a(l, n)^, 

Ji(n) = (a(0,n) + u;a(l,n))(a(0,n+l)-a;oa(0,n)). (60) 
We have decomposed the term ac as: 

ac = C2 + 2w*(T with 

6 = 2y3 . . (Jo(l) - u; • (1 - u;a:o) • Ji(l)) 

+ a;oU'^ • (Jo(2) - w {\ - wx^) ■ Ji{2)) + xgu.^ a(0, 3)^) , 

and 

00 

a^ Vl-Y^w'^xl-''. (61) 

[{wxo - 1) Ji(p + 3) - wa(l,p + 3)2 + w^xoa(0,P + 4)2y 

Collecting all previous results, and after some manipulations, x^'^^ can then be 
written as: 

X(3) = 8z.9.(xlf +2u;4.s), 

where 

xif = ^(a+C2) + ^^M + 2u;.(F(l,0) + i^(0,l)), (62) 



Series expansion method and ODE of x^^^ 19 
and 

oo 

S= <T + ^u.3p . (^i^(0,p + 2) + + 2,0)) 

p=0 

oo oo 

+ 2^^w^^+^P ■ F{k + l,p+l). (63) 

fc=0 p=0 

Appendix D gives the detailed derivation of the closed form (|27|l from H62() . 
8. Appendix C 

The purpose in this Appendix is to derive relation H58|) . From the definition of the 
variables Xn and Z„ = exp(i0„), one can deduce the following identity: 

17 7\ "^^2:1X2 Z1Z2 

(Zi - Z2) = w • (x2 - xi) — — . (64) 

1 — W^XxX^ 1 — -^1^2 

The first step is to simplify the expression of Giip) defined in (|57|) . Taking the 
integral representation of a{k,n) and v(k,n) (n = p, p + 1, p + 2), and summing on 
the index k, the sums crj(p) become: 



cro(p = u: P < Z/i?;2-; 5 -; 5 ^1^2 ^ ^1^2 - ^1 >, 

1 — w'^2:ia:2 1 — w^xiX2 



<yiy^-, — zz^- — .,2 ^ ^ iZiZ2nzi^wxo){Z2-Zi) >), (65) 



a,{p) = w-^P < v,V2 ^ ""^f ^ (ZiZ2)P(ZiZ2 - >, 

1 — w^xiX2 1 — i«a:2 

1 -2p / a:iX2 1 + ^X2, , , . 
c^2(p) = ttW p < ?;i?/2-; 2 {ZiZ2f{Zi - Z2)> . 

2 1 — W"'XlX2 1 — 10X2 

The sum C2{p) (sum of the cr's) becomes: 

C2(p) = W'^P ■ (< yiy2 ,^ ™ .2 ^Z,Z2n-\{Z2 - Zl)^) > + 

v (1 — 10^0:10:2)"' 2 

l + wa:2 X1X2 

' 1 — WX2 (1 — 10^2:1X2) ' 

where, for the first term in the brackets at the right hand side of (|65|l . the quantity 
Z1Z2 — Z\ has been replaced by —^{Z2 — Zi)^ due to the symmetry of the rest of the 
integrand in the integration variables. Using H64() and the following relation: 

<^Z"+i>^ w''+^ail,n+l), n>0 

X 

then expanding in the Z variables, working out the integration over the angles and 
making use of the identities: 

l + wx 1 + WX p 

< y cos(p</>) > = < y ZP > = yox^ujP, (66) 

1 — wx 1 — wx 

viO,p) ^ ^■{yoxl-a{0,p)), (67) 



one obtains 



C2ip) = Soip) + Slip) + E yo^o^' + P)' (68) 
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with 

Soip) - ^o(j +P + 1), Slip) = E 

Jo and Ji are defined in (|60|l . wliile J2 reads: 

J2{tT') — xq ^a(0, n) — i/7a(l, n)^ — ^a(0, n + 1) — u;a(l, n + 1)^ . 

C2{p) is now written as a sum of contributions quadratic in a's and linear in a's (the 
last term). The latter will be shown to cancel Ci{p) identically in (|56|l . Using the 
following relation: 

) -w'^a{0,n+l) (69) 

one obtains the identity: 

J2in) = ^^{Ci{n-l)-w''xoCi{n)y 
The last sum in (|58|l simply reads: 

J2w^^''yoxi+''-Mj+p) = lyo4c,ip). 

Collecting terms in H56|l . one gets: 

C{p) = yoa;g- (5o(p) + 5i(p)) (70) 

Some manipulations, such as shifts in the indices of the sums, allow us to evaluate 
two summations and to cast the quantity of interest (|58|l as: 

00 

C(0) +2^ u;3PC(p) ^ yoCiO) -2wi/o-io +10- (71) 
p=i 

^0 is a sum that can be worked out using l(T5|l . ifTTji . l(5T|l and (jS2): 

ri=0 

This ends the derivation of (|58|l given in Appendix B. 
9. Appendix D 

In the closed part Ht)2|) , there is still a summation to be done in the term ^1 defined in 
H59|l . namely the evaluation of C(0). For this purpose, we consider C{p) in the form 
GOJ). 

To evaluate 5*0(0), one shifts the index of summation to cast it as: 

^o(O) - 51 - 1) Mn) - So{0) - fo, (73) 

n=l 

with 

00 

5^(0) = nw^"-'Mn). 

n=l 

At this point, note that x" has the following Fourier expansion: 

C30 

= b{n,0) +2Ew''6(n, fc) cosfc^, 
fc=i 
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where b{k, n) reads (with m = k + n): 

bik,n)= (";^)- (74) 
„ / (1 + to) (1 + m) (2 + to) m , 

and satisfies: 

b{i,p) — a{i — l,p) — w^a{i + l,p), «>! (75) 

n 6(1, n) — b{n, 1). 

Using H75|) and the symmetry of a{k,p) in fc and p, S'o(O) becomes: 

5^(0) = 2 H ^^""^ ■ ^)^(^' 1) ~ 2Wa(n, 1)^ 

ri=l 

+n a(n, 0) • (a(n, 0) + u;^a(n, 2)^ , 
and, in integral form, we obtain: 

-50(0) = Tr^< in-, 2 cos(02) > 

ZW 1 — W X1X2 

1 ■u;^a;ia;2 „ r , \ , , : \ r , w 
+ TT^ < V1V2Y-, 2 tt2cos(02) (C0S((?!)2) - cos(0i)) > . 

Using the identity (|64ll and relation (|15|l . 5*0(0) simply becomes: 

^(0) - :^(yiy2X2Cos(</)2)) = ia(0,0)a(l,l). 
For the evaluation of S\{p\ we make use of (|69|) to get: 

S'i(0)= (a(0,^^+l)• (a(0,n-l) + w2a(0,n+l)) 

Tl=l 

+ a(0, n) ■ {2xo'w{a{0, n) + wa{l, n)) - (1 + 2w)a(0, 71 + 1)) 
Using l(T7|) . l|nU, lESl and the useful identities: 

C30 

< cos{(j)) >^ 2w a(0, 0)a(0, 1) + 2 ^ a(0, n) a(0, n + 1), 

00 

< y2 cos(20) > = w^a{0, 1)^ + 2 ^ w^" a(0, n - 1) a(0, n + 1), 



we obtain 



n=l 



SliO) = TT- ((a(0, 0) + wa(l, 0)) • a(0, 0) - < (1 + li;^;) >) 

ZW 

+ < y'^ cos((?!)) (1 - 2w cos{(j)) + 2w) > 

+ ^(3^^0(0, if + a(0, of -2(1 + 2w) a(0, 0) a(0, 1)). 



With the relation p5|l and using Ht)7|l . we obtain 
< 2/^ • (1 + wx) > = -y(0, 0) = 



2w 
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With the following relation: 

1 — 2 w cos(0) = h w'^x, 

X 

and using (|66|l . we get: 

< cos(0) (1 — 2'u;cos(0) + 2w) > = wxiqXq 

and finally: 

S'i(0)= ^(xoa(0,0) + (l + 2wa;o)a(0,0)2 (76) 

+ 3^^0(0, 1)^ + 2{w^X(i - 1 - 2w) • a(0, 0)a(0, 1)) . 

The last step in calculating the closed form (|62|l is to explicitly evaluate the terms 
F(0, 0), 0) and i^(0, 1). This can be done with the help of Appendix A. Note that 
g{\, 1), and g{0, 0), can be written as: 

5(1, 1) = ^ (5(1' 0) - 2g(0, 0) + 2v(0, 0) - a(l, 0)) , 

E and K being defined in ^F!^ . 

Collecting all the previous results, expanding in the basis of E and K functions, 
one obtains the relation (|27|l . where the polynomials qij and pij are given by: 

(?oo = - 96u;^ + 276u;^ - 120w'^ + 144^^ - 64u;2 + 2,2w - 16, 

gio = 128w^ - 832zx;'^ - 792m;^ + 740u;'' - 680^^ + 295^^ 

+ 32w- 16, 

(701 = 68w^ - 212u;^ - 120w;^ + 89w^ - 9Qw + 48, 
(720 = 960'u;^ - 416u;^ + 1924w5 - 1048w'^ - 170^^ + 98u;2 
+ 3u; - 2, 

(711 584w5 - 256u;'* + 7Uw^ - 386u;2 - 38u; + 20, 

(702 = 82^3 - i2v? + 67u; - 34, 

530 = {Aw - if {Aw + l)2(2w - l)(2w + l)w^ 

(721 = - {Aw- l){Aw + l){Aw^ - 28w^ - 2w + 5)w^, 

(712 = - 8 (6w^ -9w^ -w + l)w^ 

(/03 - -2(3u;-2)«;2, 



Poo = 528w'^ - A6Aw^ + 2A8w^ - 212w^ + 96u;2 - 64u; + 16, 
Pio = - 1280u;^ + 352w'^ + 1984u;5 - 1504u;'* + 1284^;^ - 264^2 

- 64w + 16, 

poi = - 320u;^ + 208w'' + 284^^ - \Mw^ + 192u; - 48, 
P20 = 256w* - 1280w^ + 2320u;^ - 4144w^ + 928^-^ + 372u;3 

- 98u;^ - 7u; + 2, 

pii = 384w6 - 896u;5 + 960?i;'^ - 1548u;3 + 356u;2 + 78w - 20, 
P02 = Sew'^ - 120w^ + 94^2 - 135w + 34. 
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10. Appendix E 

This Appendix contains the details of evaluation of x^^^ using the yx" expansion 
method. One should note however, that x'^^ can be obtained in a straightforward 
manner by direct integration. 
One writes x^^'' as (7|: 

X^'^ = (77) 



where 

x(^).4«.^.<,V^i±^.(l-cos(20))>. 

One defines: 

A{k,n) = A{-k,n) = w^''^ a{\k\ , n) 
With the variable Z — exp(i0), the expansion of yx" is written as: 

oo oo 

yx" ^ J2 A{k,n)Z^ = ^ A{k,n)Z-^ 

k — — OQ k — — OQ 

Expanding in a;, x^'^^ becomes : 

oo 

X(2) = Aw^ ^(n + 1)2 it,2" < y22.2n+2(^ _ cos(20)) > . (78) 
n=0 

Now write < y'^x^^+'^ (1 - cos(2(/))) > as < ya;"^ (1 - cos(2(/))) > with 

ni + n2 = 2n + 2. One has: 

< - cos(2(/))) > (79) 

OO 

= A(fci,ni)A(fc2,n2) < (1-^2/2 -Z-2/2)Z^-i+'=^ > 

ki ,k2— — oo 

OO 

= J2 ^(fc,«i)(^(fc,"2)- A(fc + 2,n2)/2- A(fc-2,?i2)/2). 

fc— — oo 

Some manipulations give: 

< yx"! • (1 - cos(20)) > (80) 

OO 

= (A(fc - 1, m) - A(fc + 1, m)) (A(fc - 1, na) - A(fc + 1, n2)) . 
fe=i 

Coming back to the definition (|78|l of A{k,n), one recognizes in (I8U|I the function 
h{k,n) defined in ifT^ . x*^^\ then, becomes (recall that ni + ^2 = 2n + 2) 

C30 oo 

x(2) = 4^/;^.^^(^^ + l)2^i;2"+2fe5(fc + l^„^)5(fc + l^„2). 

One thus has many equivalent forms for x*^^^ depending on the partition ni +ri2 = 
2n + 2. The above sums merges into a single hypergeometric function: 

X''^ = 4^4.2Fi(y;3;16«;^), (81) 
which is the well-known result 1 : 
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Here we give an alternative derivation of (|81(l . performed, this time, with another 
expansion of the quantity yx", which underhnes the role played by some fusion 
relations on hypergeometric functions, 
j/x" can be written as : 



yx" 



3=0 



n + j 
n 



sm(-) ^ • 



2Fi{n+l/2,n + j + 1; 2n + 1; Aw 



(83) 



Taking x^^^ in the form (|78|) . writing y^x^""''^ as ya;"+^ ya;""''^, using the 
expansion H83|) for each yx", the angular integration gives: 



Aw" 



3 

E C(ri, fci,fc2) • 2i^i(n + -,n + 2 + fci;2n + 3;4w) 



4w 



(84) 



where 



C(n, fci, fc2) = U;2n+fci+fc2 . (_4)fcl+fe2 . („ + 1)2 (3/2)fei+fc2 



(3)fe, 



+ fe2 



1 + fcl 



n + l + k2 
k2 



and (-/V)„ are the usual Pochhammer symbol. It is worth comparing H81(l with l|84(l . 
One sees, with this alternative expression for x^^^ that the complexity of x*-^' 
now been encapsulated in a single hypergeometric expression 2^"i (^5/2, 3/2; 3; 16 
In the holonomic approach developed in this paper, the occurrence of the prototype 
of holonomic functions, namely hypergeometric functions, is quite natural. These 
calculations also underline the important role the fusion identities |46) (see relation 
(1.9) in '46') on hypergeometric functions, like F = ■ F F, can play in such 

calculations. This is probably the sign of quite deep structures, and symmetries, like 
the fusion relations encountered in CFT, or integrable models, and deeply connected 
with Yang-Baxter relations "47", '48' . 

At this point it is worth recalling that the fusion relation (|84ll is totally reminiscent 
of, for instance, the fusion relation (1.9) in ' 



Or, Ct 



bs, du 



r+tFs-\-u 

X k+tFm+u+1 



; xw 



E 



(Q)n(efc)„(~a;)" 
{du)n{fm)n{n + 7)„ri! 



m+r+2 



Fk+s 



n=0 

n + Ct, n + ek 
2n + 1 + 7, n + du, n + f, 
-n,n + -f, Or, /, 

bs, ek 



w 



where 7, e and / are arbitrary and ar, ct, ■ ■ ■ are r, t, ■ ■ ■ upper and lower parameters. 

This fusion relation, of the type F{uv) = ^ C ■ F{u) F{v), gives a representation 
of a hypergeometric function of the product uv, as a sum of products of hypergeometric 
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functions of u and hypergeometric functions of u. In this respect relation (|84|l 
corresponds to particular cases of fusion relations F{uv) — ' with 

u = V = Aw. 



12. Appendix G 

To evaluate the sums in H25|) in order to obtain our long series for x''^\ substantial 
computing time is spent in the computation of the coefficients of the various functions 
a{k,p), g{k,p), v{k,p), and i/qXq. While the series expansion of a(k,p), and i/oXq, 
simply read 



oo 



yoxl 



■_ V j J \ j + k 

j=0 \ J / 

the coefficients of the auxiliary function g{k,p) are given by recurrences of depth three 
with huge polynomials in fc, p and j. For series generation purposes, we have found it 
more efficient to use the following inhomogeneous recurrences between the coefficients 
of g{k,n) and a{k,n). 

Noting by Cg{k,p,i) the coefficient of in the expansion of g(fc,p) (and similarly 
for a{k,p) and v{k,p)), the inhomogeneous recurrences read: 

0= Cg{n + 2,p,i) -2Cg{n + l,p,i + l) + Cg{n,p,i + 2)- 

Cain+l,p,i/2+l), I even , , 

0, iodd. ^^^> 

Once the coefficients of g{k,p) are found, those of v{k,p) are given by 

Cy{n,p,i) = Cg{n,p,i)-Cg{n+l,p,i-'l). 

The following identities are useful for (|87|l : 

Cgik,p,0)^ Ca{k + l,p,0), Cg{k,p,l) = 2Ca{k + 2,p,0), (88) 

Cgil,p,^)=CgiO,p,^+l) -K^^'^l^^^ 

i even 

(89) 

iodd 

where is deduced from ^7^. 




13. Appendix H 

Recall that the separation of x^^-* into a relatively simple closed form H26|l . a "hard 
to compute" series S/i, and a simpler one Ss, is made for algorithmic considerations. 
Although this separation is not " natural" , it is tempting to seek immediately, as 
a first step, a linear differential equation for (|26|l which should be "simpler" than 
X*^^'. Having two functions, / and g, satisfying two homogeneous linear differential 
equations (like the second order differential equations satisfied by K and E (see 
((23), there are some formal computer programs (gfun or others 49 ) which enable 
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one to build the homogeneous linear differential equation satisfied by the product 
f ■ g or the sum f + g, and of course step-by-step, a (cubic) sum of products like 
H26|l . For instance, from the two order- two differential equations for K and E, one 
builds the order four differential equation of the product K ■ E and the order six 
differential equation of the cubic term K ■ E'^ . The latter term (with a polynomial 
multiplying the highest derivative of degree 17), combined with its associated algebraic 
expression, gives an order six differential equation with a polynomial multiplying the 
highest derivative of degree 38. This step by step procedure gives more and more 
complicated linear differential equations. It seems pointless to give these extremely 
involved, and cumbersome, expressions. Just note that the linear differential equations 
corresponding to the sum of the strictly cubic terms in E and K will be of order 20, 
that the one for the strictly quadratic terms in E and K will be of order 18, of order 
8 for the linear terms in E and and, of course, of order one for the algebraic part 
in H26|) . yielding a linear differential equation for (|26|) of order 47 and an extremely 
large degree 267 as a lower bound| This means that using a bruteforce approach based 
on series expansions, instead of the use of the closed form, would have required more 
than 12500 coefficients for finding this homogeneous linear differential equation! As 
a consequence, it might look hopeless, at first sight, to seek a homogeneous linear 
differential equation for x^"^"*- 

From the series point of view, let us focus on the growth of the coeflacients of the 
three series comprising x^'^^- The integer coefficients of w'^^^ read, respectively, for 
the "closed", "simple", "hard" and x^'^^ /divj^ series : 

Cci ^ 24752.108 •4'^^\ Cs ~ -25012.475 • 4''^\ 

Ch ^ 271.423- 4'**^\ C ~ 11.056 • 4^^^ 

One thus sees that the coefficients for x^^^ ^-re actually obtained from mutual 
cancellations of much larger integers: the coefficients of Xc/ ^ ^^^d (w) , are of the 
same order, but of opposite signs, their sum, Cd + Cg — —260.366 • 4"'®^, is of the 
same order as the coefficient of 'E,fi{w), but of opposite sign, the total sum becoming 
the coefficient of x''^''- The coefficient of w^^^ for the closed form Xc?'' is — 2238.72 
times larger than the corresponding coefficient for x^"^^! 

These points show that x^^^ seen as the sum (I25II . should be much simpler than 
the sum of its constituents. This is indeed what was discovered. 

[1] T.T. Wu, B.M. McCoy, C.A. Tracy and E. Barouch, The spin-spin correlation function for 
the two-dimensional Ising model: exact results in the scaling region, 1976 Phys. Rev. B 13, 
316-374 

[2] C. R. Nappi, 1978 Nuovo Cim. A 44 392 

[3] J. Palmer, C. Tracy, 1981 Adv. Appl. Math. 2 329 

[4] K. Yamada, 1984 Prog. Theor. Phys. 71 1416 

[5] K. Yamada, 1985 Phys. Lett. A 112 456-458 

[6] B. Nickel, 1999 J. Phys. A: Math. Gen. 32 3889-3906 

[7] B. Nickel, 2000 J. Phys. A: Math. Gen. 33 1693-1711 

[8] N. Zenine, S. Boukraa, S. Hassani, J.M. Maillard, The Fuchsian differential equation of 
the square lattice Ising model x'^' susceptibility, 2004 J. Phys. A 37 9651-9668 and 
arXiv: math-ph /0407060 
[9] I. Syozi, S. Naya, 1960 Prog. Theor. Phys. 24 829 
[10] D. Hansel and J-M. Maillard, 1988 J. Phys. A: Math. Gen. 21 213-225 
[11] D. Dhar and J-M. Maillard, 1985 J. Phys. A: Math. Gen. 18 L383-L388 



f Much larger than in the seventh order Fuchsian equation given in (8j for x'^' and recalled in the 
following (sec I34i above). 



Series expansion method and ODE of x^^^ 



27 



[12] R.J. Baxter, 1982 J. Stat. Phys. 28, 1-41 

[13] M.T. Jaekel and J-M. Maillard, 1985 J. Phys. A: Math. Gen. 18 641, 1229, 2271 

[14] D. Hansel and J-M. Maillard, Formal constraints on series analysis on the Potts models, 1987 

Mod. Phys. Lett. B 1, 145-153 
[15] D. Hansel and J-M. Maillard, Series analysis on the q-state checkerboard Potts models, 1988 

Int. Journ. Mod. Phys. B 1447-1462 
[16] C. N. Yang, Phys. Rev. 1952 85, 808-816 
[17] C. H. Chang, Phys. Rev. 1952 88, 1422 

[18] D. Hansel, J.M. Maillard, J. Oitmaa, M.J. Vergakis, J. Stat. Phys. 1987 48 69-80 

[19] A.J. Guttman, I.G. Enting, 1996 Phys. Rev. Lett. 76 344-347 

[20] L. Lipshitz, 1989 J. Algebra 122 354 

[21] D. Zeilberger, 1990 J. Comp. Appl. Math. 32 321 

[22] P. Cartier, Fonctions polylogarithmes, nombres polyzetas et groupes pro-unipotents, Asterisque 

282 (2002), 137-173 (Sem. Bourbaki no. 885). 
[23] T. Terasoma, Mixed Tate motives and multiple zeta values, 2002 Invent. Math. 149, 339-369; 

preprint AG/0104231. 

[24] M. E. Hoffman and Y. Ohno, Relations of multiple zeta values and their algebraic expression, 

2003 J. Algebra 262, 332-347; preprint QA/0010140. 
[25] S. Oi, Representation of the Gauss hypergeometric function by multiple poly logarithms and 

relations of multiple zeta values, preprint NT/0405162 
[26] H.K. Ince, Ordinary differential equations, (Longmans, London, 1927; also Dover Pubs., NY, 

1956) 

[27] U. Mugan and A. Sakka, Schlesinger transformations for Painleve VI equation, 1995 J. Math. 
Phys. 36, 1284-1298. 

[28] P. Painleve, Memoire sur les equations differentielles dont I'integrale generale est uniforme, 

1900 Bull. Soc. Math. Phys. France 28, 201-261. 
[29] A. V. Kitaev, Special Functions of the Isomonodromy Type, Rational Transformations of 

Spectral Parameter, and Algebraic Solutions of the Sixth Painleve Equation, 

http://arXiv.org/abs/nlin/0102020 
[30] E. Barouch, B.M. McCoy and T.T. Wu, 1973 Phys. Rev. Lett. 31, 1409-1411 
[31] J. Chazy, Sur les equations differentielles dont I' integrale possede une coupure essentielle 

mobile, 1910 C. R. Acad. Sc. Paris 150, 456-458 
[32] J. Chazy, Sur les equations differentielles du troisieme et d' ordre superieur dont V integrale 

generale a ses points critiques fixes., 1911 Acta Math. 34, 317-385 
[33] R. Conte, The Painleve Property, One Century Later, CRM Series in Mathematical Physics, 

Springer, R. Conte editor. 
[34] M. D. Kruskal, N. Joshi and R. Halburd, Analytic and asymptotic methods for nonlinear 

singularity analysis: a review and extensions of tests for the Painleve property, (1996) 
[35] W.P. Orrick, B.C. Nickel, A.J. Guttmann, J.H.H. Perk, 2001 J. Stat. Phys. 102 795-841 
[36] W.P. Orrick, B.C. Nickel, A.J. Guttmann, J.H.H. Perk, 2001 Phys. Rev. Lett. 86 4120-4123 
[37] B.M. McCoy, T.T. Wu, 1980 Phys. Rev. Lett. 45 675-678 
[38] J.H.H. Perk, 1980 Phys. Lett. A79 3-5 

[39] M. Jimbo, T. Miwa, 1980 Proc. Japan Acad. A 56 405; 1981 Proc. Japan Acad. A 57 347 
[40] J J Rehr, G S Joyce and A J Guttmann, A recurrence technique for confluent singularity analysis 

of power series, 1980 J. Phys. A: Math. Gen. 13, 1587-1602 
[41] M. H. Khan, High-order differential approximants, 2002 Journal of Computational and Applied 

Mathematics, 149, Issue 2, 457-468 
[42] T. Sasaki and M. Yoshida, A system of differential equations in 4 variables of rank 5 invariant 

under the Weyl group of type Eq, 2000 Kobe J. Math. 17, 29-57. 
[43] J. M. Maillard, Hyperbolic Coxeter groups, symmetry group invariants for lattice models in 

statistical mechanics and Tutte-Beraha numbers. 1995 Math. Comput. Modelling 26, 169- 

225 

[44] S. Boukraa, J-M. Maillard, and G. RoUet, Discrete Symmetry Groups of Vertex Models in 

Statistical Mechanics. 1995 J.Stat.Phys. 78, 1195-1251. 
[45] S. Boukraa, and J-M. Maillard, Let's Baxterise. 2001 J. Stat. Phys. 102, 641-700 and : 

hcp-th/0003212 

[46] G. Gasper, Using symbolic computer algebraic systems to derive formulas involving orthogonal 
polynomials and other special functions, published in Orthogonal Polynomials: Theory and 
Practice, ed. by P. Nevai, Kluwer Academic Publishers, Boston (1990), pp. 163-179. 

[47] J. Avan, J-M. Maillard, M. Talon and C-M. Viallct, New local relations for lattice models. 1990 
Int. Journ. Mod. Phys. B 4, 1895-1912 



Series expansion method and ODE of x^^^ 



28 



[48] C. Gomez, M. Ruiz-Altaba, G. Sierra, Chapter 9 Duality in conformal field theories in Quantum 
Groups in Two-dimensional Physics, Cambridge monographs on Mathematical Physics, 
Cambridge university press, 1996 

[49] Mgfun's project : see http:/ /algo.inria.fr/chyzak/ mgfun.html| ; gfun - generating functions 
package see gfun in : ,http: / / algo.inria.fr /libraries! 



